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Abstract 

We measure the neutral D total forward cross section and the differential cross 
sections as functions of Feynman-a; (xf) and transverse momentum squared for 
500 GeV/c 7r~-nucleon interactions. The results are obtained from 88 990±460 
reconstructed neutral D mesons from Fermilab experiment E791 using the decay 
channels D -^K~n + and D — > K^tt + tt^tt + (and charge conjugates). We extract 
fit parameters from the differential cross sections and provide the first direct mea- 
surement of the turnover point in the xf distribution, 0.0131±0.0038. We measure 
an absolute D +D (xf>0) cross section of 15.4 i 2 J /ibarns/nucleon (assuming a 
linear A dependence). The differential and total forward cross sections are compared 
to theoretical predictions and to results of previous experiments. 
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Charm hadroproduction is a convolution of short range processes that can be 
calculated in perturbative quantum chromodynamics (QCD) and long range 
processes that cannot be treated perturbatively and thus must be modeled us- 
ing experimental measurements. The large theoretical uncertainties from both 
contributions are reflected in the relatively large number of input parameters 
that can be adjusted when comparing models to the results of experiments. A 
single measurement, no matter how precise, cannot unambiguously determine 
these parameters. However, the results of high statistics measurements like the 
ones reported here, when combined with other measurements of similar pre- 
cision, can constrain such parameters as the charm quark mass, the intrinsic 
transverse momentum of the partons in the incoming hadrons, and the effec- 



tive factorization and renormalization scales used in theoretical calculations. 

We report here measurements of the differential cross sections versus the kine- 
matic variables Feynman-x (xf) and transverse momentum squared (p^), as 
well as the total forward cross section for the hadroproduction of neutral D 
mesons. The relatively high pion beam momentum, 500 GeV/c, coupled with 
the good geometric acceptance of the Tagged Photon Laboratory (TPL) spec- 
trometer, allows us to investigate a wide kinematic region that includes points 
at negative xp. We are able to measure the shape of the differential cross 
section versus xp with sufficient precision to confirm, for the first time, that 
the turnover in the cross section does occur at xf>0, as expected for incident 
pions [1]. 

Combining data from two D decay modes, D — > Kn and D — ► Knim, Q 
we extract a sample of 88 990±460 (78730±430 at xf>0) fully reconstructed 
charm decays to use for these measurements. In addition to the greater statis- 
tical significance, the use of two modes provides a means to better understand 
the systematic errors associated with the reconstruction of the decay products 
of these fully-charged decays. 

The data were accumulated during the 1991/1992 Fermilab fixed-target run of 
experiment E791 [2,3]. The experiment utilized the spectrometer built by the 
previous TPL experiments, E516 [4], E691 [5], and E769 [6], with significant 
improvements. The experiment employed a 500 GeV/c tt~ beam tracked by 
eight planes of proportional wire chambers (PWC's) and six planes of silicon 
microstrip detectors (SMD's). The beam impinged on one 0.52-mm thick plat- 
inum foil (1.6 cm in diameter) followed by four 1.56-mm thick diamond foils 
(1.4 cm in diameter), each foil center separated from the next by an average of 
1.53 cm, allowing most charm particles to decay in air. The downstream spec- 
trometer consisted of 17 planes of SMD's for vertexing and tracking along with 
35 planes of drift chambers, 2 PWC planes, and 2 analysis magnets (bending 
in the same direction) for track and momentum measurement. Two multi- 
cell threshold Cerenkov counters, an electromagnetic calorimeter, a hadronic 
calorimeter, and a wall of scintillation counters for muon detection provided 
particle identification. The trigger was generated using signals from scintilla- 
tion counters as well as the electromagnetic and hadronic calorimeters. The 
beam scintillation counters included a beam counter 1.3 cm in diameter (14 cm 
upstream of the first target) and a large beam-halo veto counter with a 1.0 cm 
hole (8 cm upstream of the first target). The interaction counter was located 
2.0 cm downstream of the last target and 0.6 cm upstream of the first SMD 
plane. The first-level trigger required a signal corresponding to at least 1/2 of 
that expected for a minimum ionizing particle (MIP) in the beam counter, no 



1 Charge conjugates are always implied. We use D to represent the sum of D and 

D . Similarly, Kir (Kinnr) includes K~tt + (K~ir + ir~ir + ) and charge conjugate. 



signal greater than 1/2 of a MIP in the beam halo counter, and a signal cor- 
responding to greater than ~4.5 MIP's in the interaction counter (consistent 
with a hadronic interaction in one of the targets). The second-level trigger re- 
quired more than 3 GeV of transverse energy in the calorimeters. Additional 
requirements eliminated events with multiple beam particles. A fast data ac- 
quisition system [7] collected data at rates up to 30 Mbyte/s with 50 /xs/event 
deadtime. Over 2xl0 10 events were written to 24 000 8mm magnetic tapes 
during a six-month period. 

The raw data were reconstructed and filtered [3,8] to keep events with at least 
two separated vertices, consistent with a primary interaction and a charm 
particle decay. Following the event reconstruction and filtering, selection cri- 
teria for the D candidates were determined by maximizing S/y/S + B where 
S is the (normalized) number of signal events resulting from a Monte Carlo 
simulation and B is the number of background events appearing in the data 
sidebands of the reconstructed Kw or Ktttttt mass distribution. Only selection 
variables that are well modeled by the Monte Carlo simulation were used. 
The final selection criteria varied by decay type (K7r and Kiririr) and by xp 
region. The full range of a cut variation is given in the descriptions below. 
To eliminate generic hadronic interaction backgrounds as well as secondary 
interactions, the secondary vertex was required to be longitudinally separated 
from the primary vertex by more than 8-11 times the measurement uncer- 
tainty on the longitudinal separation (~400 /im) and to lie outside of the 
target foils. Backgrounds from the primary interaction were also reduced by 
requiring that the candidate decay tracks miss the primary vertex by at least 
20-40 /im. To ensure a correctly reconstructed charm particle and primary 
vertex, the momentum vector of the D candidate was required to point back 
to within 35-60 /im of the primary vertex and to have a momentum compo- 
nent perpendicular to the line connecting the primary and secondary vertices 
of less than 350-450 MeV/c. Finally, the sum of the squares of the transverse 
momenta of the decay tracks relative to the candidate D momentum vector 
was required to be greater than 0.4 (GeV/c) 2 (0.15 (GeV/c) 2 ) for the Kit 
(Kinnr) candidates to favor the decay of a high-mass particle. All primary 
vertices were required to occur in the diamond targets. Thus, our results come 
from a light, isoscalar target. The Cerenkov information is not used in this 
analysis; all particle-identification combinations are tried. The inclusive Kit 
and Kkiht signals are shown in Fig. 1. 

The reconstructed data were split into 20 bins of xp, integrating over all p T , 
and 20 bins of p^, integrating over xp>0. To combine data from varied con- 
ditions (e.g., using particles that pass through one and two magnets), the 
normalized mass (m n ) is constructed for each candidate using its calculated 
mass and error (m and a m ) and the measured mean mass (m^): m n = m-mg . 

v / \ / (J m 

Using the binned maximum likelihood method, the normalized mass distri- 
butions were fit to a simple Gaussian for the signal and linear or quadratic 
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Fig. 1. The Kir and Kiririr signals for data with — 0.125<Xi?<0.8. 
polynomials for the background. 

The acceptance can be factorized into trigger efficiency (e tr ig) and reconstruc- 
tion efficiency (e rec ). Most of the trigger inefficiency is due to vetoes on multiple 
beam particles. These resulted in a (70.3±1.1±4.2)% trigger efficiency, where 
the first error is statistical and the second error is systematic. The interaction 
and transverse energy requirements were greater than 99% efficient for recon- 
structable hadronic charm decays. Writing and reading the data tapes was 
(97.6±1.0)% efficient. Combining these efficiencies gives e trig = (68.3 ±4.4)%, 
where the error is dominated by the systematic error. The reconstruction 
efficiency is obtained from a Monte Carlo simulation. The Monte Carlo simu- 
lation used Pythia/Jetset [9] as a physics generator and models the effects 
of resolution, geometry, magnetic fields, and detector efficiencies as well as all 
analysis cuts. The efficiencies were separately modeled for five evenly spaced 
temporal periods during the experiment. This was motivated by a highly in- 
efficient region of slowly increasing size in the center of the drift chambers 
caused by the 2 MHz pion beam. The Monte Carlo events were weighted to 
match the observed data distributions of xp, p%, and the summed p\ of all 
two-magnet charged tracks in the event other than those from the candidate 
D meson. The resulting reconstruction efficiencies as a function of xf and p\ 
are shown in Fig. 2. 

From the number of reconstructed D candidates, the reconstruction effi- 
ciency, the trigger efficiency, and the PDG branching fractions [10], we obtain 
the number of D mesons produced in our experiment during the experiment 
livetime, N pro d- The cross section as a function of each variable z (where z = 
Xp or j?y), is: 



a(ir-N -► D v X:z) 



N prod (D°;z) 

T N N n - 



Tjy is the number of nucleons per area in the target (calculated from the target 
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Fig. 2. D — ► ivr-7r and D — ► Ktx-ktx reconstruction efficiencies versus xf (left) and 
p T (right). 

thickness) and is (1.224 ± 0.004) x 10~ 6 nucleons//xb. N w - is the number of 
incident 7r~ particles during the experiment livetime. This is obtained directly 
from a scaler which counted clean beam particles (>l/2 MIP signals in the 
beam and interaction counters and no signal greater than 1/2 MIP in the 
beam-halo veto counter) during the experiment livetime. These are the only 
beam particles which could cause a first-level trigger. Using Eq. 1, we obtained 
the D +D differential cross sections versus xf and p T shown in Figs. 3 and 
4. 



The systematic errors are divided into two categories and incorporated in two 
stages. The uncorrelated systematic errors are determined individually for the 
Kit and K7nnr results. These systematic errors include uncertainties in the 
Monte Carlo modeling of the selection criteria, the background functions, and 
the widths used in the Gaussian signal functions. The correlated systematic 
errors are calculated for the combined D result, obtained from adding the 
Ktt and Kinrir samples together, weighted by the inverse- square of the com- 
bined statistical and uncorrelated systematic errors. The correlated errors are 
associated with uncertainties in the D lifetime, the Monte Carlo production 
model, the Monte Carlo weighting procedure, and the run period weighting 
procedure. Finally, we compare our measured Kn to K-rnnr branching ratio to 
the PDG [10] value to estimate the residual tracking and vertexing efficiency 
modeling error. In addition to these errors, which can affect both the normal- 
ization and the shape of the differential cross sections, there are two errors 
which affect only the normalization: the uncertainties in the trigger efficiency 
and target thickness. 

For the differential cross sections shown in Figs. 3 and 4, the systematic errors 
are factorized into shape and normalization parts. The error bars in the figures 
show the sum, in quadrature, of the statistical and all systematic errors after 
factoring out the normalization component. Although the relative importance 
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Fig. 3. Fits to the -D +.D xp differential cross section with the functions given in 
Eqs. 2 (dashed) and 3 (dotted). Error bars do not include a i"i§% normalization 
uncertainty. 



varies bin-by-bin, the most important systematic errors generally come from 
uncertainties in the signal width and the Monte Carlo efficiency modeling. In 
all cases, the systematic error dominates. For the total forward cross section, 
all of the errors are summarized and summed in Table 1. 

In the past, xp distributions have been fit with 



da 

dxp 



A(l-\x F \) n . 



(2) 



Fitting Eq. 2 in the range 0.05 <xp< 0.50, we find n = 4.61±0.19, as shown in 
Fig. 3. This function does not provide a complete representation of our data. 
Although the x 2 /dof is small (0.3), the value of n is quite dependent on the 
range fitted and on the errors on the data points. Another function, which can 
be extended into the negative xp region, is an extension of Eq. 2 which uses 
a power-law function in the tail region and a Gaussian in the central region; 



+. 



o 






-a 



*■* 



•"^ 



• E791 D data (K71 + Kjt7tJt) 

Fit to exp(-bp-f-) (0<p|<4) x 2 /dof=2.8 

b = 0.83 ± 0.02 

Fit to (.ami + Pt) ~ P (0<p|<18) x 2 /dof=0.3 

(X = 2.36 ±0.23 P = 5.94 ±0.39 



Pt 



3.5 

(GeV/c) 2 



• E791 D data (K71 + KjKtJl) 

Fit to exp(-b 'p T ) ( 1 <pr< 18) 

b' = 2.41 ±0.03 



o 



2.36 ±0.23 



(0<pf<18) 
P = 5.94 ±0.39 



X /dof = 1.7 
X 2 /dof =0.3 



2 

Pt 



16 18 

(GeV/c) 2 



>0 , t=:0 



Fig. 4. Fits to the D u +D u p\ differential cross section with the functions given in 
Eqs. 4 (dashed, top), 5 (dashed, bottom), and 6 (dotted, top and bottom). The top 
plot shows the range 0<p^<4 (GeV/c) 2 while the bottom plot shows the full range, 
0<p|,<18 (GeV/c) 2 . Error bars do not include a ±"15% normalization uncertainty. 



that is, 



da 

dXF 




x F -x c \) n ', 

1 / x F -x c \2 
2\ cr I 



\Xf 

\xf 



x c 



> x b 
< x b 



(3) 



Requiring continuous functions and derivatives allows us to write Eq. 3 with 
one normalization parameter and three shape parameters: n' gives the shape 
in the tail region, x c is the turnover point, and x b is the boundary between 
the Gaussian and power-law function. The fit parameters from this func- 
tion are nearly independent of the fit range. Fitting our data in the range 
-0.125<:r F <0.50 gives n' = 4.68 ± 0.21, x c = 0.0131 ± 0.0038, and x b = 
0.062 ± 0.013 with a x 2 /dof '=0.4, as shown in Fig. 3. This is the first mea- 
surement of the turnover point x c in the charm sector. The fact that it is 
significantly greater than zero is consistent with a harder gluon distribution 
in the beam pions than in the target nucleons. 



Table 1 

Sources and values of the uncertainties on the total forward D +D cross section 
measurement. The total systematic error comes from the quadratic sum of the D 
systematic errors. Although some of the systematic errors are subdivided, each D 
systematic error is obtained directly (including correlations) and is not the quadratic 
sum of the components. 



Error type 




Error (%) 






Kir 


K-rnrir 


D° 


Statistical & uncorrelated systematic errors 


+4.7 
-9.0 


+7.3 
-11.0 


+4.7 
-7.6 


Statistics 


±1.0 


±2.8 




Selection criteria efficiency modeling 


+0.1 
-0.0 


+1.2 
-0.0 




MC background function 


+0.0 
-5.5 


+0.0 
-6.0 




D signal width in fits 


+4.0 
-4.1 


+3.9 

-4.4 




PDG [10] branching ratio 


±2.3 


±5.3 




Correlated systematic errors 






+6.9 
-9.9 


D lifetime 


+0.6 
-0.5 




Monte Carlo kinematic weighting 


+0.0 
-2.3 




Monte Carlo production model 


+5.9 
-5.1 




Time dependent efficiency modeling 


+1.0 
-2.0 




Tracking and vertex finding 






±4.6 


Trigger efficiency 






±6.4 


Target thickness 






±0.3 


Total 






+11.5 
-14.8 



The functions which have been used in the past to fit the p\ distribution are: 

da 



dp\ 



Ae -bfi 



(4) 



at low p\ {p\ <4.0 (GeV/c) 2 for this analysis), 



da 
dp\ 



Ae~ b ' PT 



(5) 



at high p\ (p\ > 1.0 (GeV/c) 2 for this analysis), and 



da 
dp\ 



A 



ami + Pt 



(6) 



over all p\ with m c set to 1.5 GeV/c 2 [11]. The results of fitting these equations 



to the data are shown in Fig. 4. For the ranges given above, the fit results are: 

• b = 0.83±0.02 with x 2 /dof = 2.8, 

• b' = 2.41±0.03 with x 2 /dof = 1.7, and 

• a = 2.36±0.23 (GeV/c 2 )- 2 and (3 = 5.94±0.39 with x 2 /dof = 0.3. 

Equation 4 does not provide a good fit even over the very limited range to 
which it is applied. While the x 2 /dof (1.7) of the fit to Eq. 5 is not good, 
it appears to be a reasonable fit to the data. Equation 6 provides a very 
good fit to the data over the entire range of p^. Unfortunately, using two 
free parameters (in addition to the normalization) makes it more difficult to 
compare to other experiments and theory since the parameters in this fit are 
highly correlated. This is reflected in the large (7-10%) errors on a and (3 
compared to the error on b' (1%), as shown above. 

Figures 5 and 6 show a comparison of our xf and p\ distributions to theoret- 
ical predictions for charm quark and D meson production. Although the data 
come from D mesons, the theoretical predictions for charm quark production 
are included for completeness. The theoretical curves are normalized to ob- 
tain the best fit (lowest x 2 /dof) to our data. The theoretical predictions come 
from a next-to-leading order (NLO) calculation by Mangano, Nason, and Ri- 
dolfi (MNR) [12] and the Pythia/ Jetset [9] event generator. The MNR NLO 
charm quark calculation uses SMRS2 [13] (HMRSB [14]) NLO parton distribu- 
tion functions for the pion (nucleon), a charm quark mass of 1.5 GeV/c 2 , and 

an average intrinsic transverse momentum of the incoming partons (W (A: 2 )) 



of 1.0 GeV/c. The value for W (/cf) was suggested by M. L. Mangano [15] and 
is independently motivated by the study of azimuthal angle correlations be- 
tween two charm particles in the same event [11]. The D meson results are 
obtained by convoluting the charm quark results with the Peterson fragmen- 
tation function [16] with e = 0.01. The low value for e was also suggested by 
M. L. Mangano [15] in response to a reanalysis of D fragmentation in e + e~ col- 
lisions [17]. The Pythia/ Jetset event generator uses leading order D02 [18] 
(CTEQ2L [19]) parton distribution functions for the pion (nucleon), a charm 



quark mass of 1.35 GeV/c 2 , J (/c 2 ) of 0.44 GeV/c, and the Lund string frag- 
mentation scheme to obtain D results. Tables 2 and 3 show a comparison of 
our xp and p\ fit results to theoretical predictions and to recent high-statistics 
charm experiments which used pion beams. The evident energy dependence 
of the shape parameters in Tables 2 and 3 are consistent with theoretical 
predictions [11]. 

We obtain the total forward cross section by summing the xp differential cross 
section for xp>0 and assuming the cross section for 0.8<a;^<1.0 is half that of 
the cross section for 0.6<xp <0.8 but with the same error. Assuming a linear 
dependence on the atomic number [22] , we obtain the neutral D total forward 
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Fig. 5. The D +D xp differential cross section compared to various theoretical 
predictions described in the text. The curves are normalized to obtain the best fit 
to the data in each case. Error bars do not include a ii?% normalization uncertainty. 

cross section, a(D + D ; x F > 0) = 15.4 l^l Atbarns/nucleon. To obtain 
the total charm cross section, a(cc), we multiply our D +D cross section by 
1.7. This accounts for three multiplicative effects: the relative production of 
charm quarks compared to D mesons (2.1), the conversion from xp>0 to all 
xp (1.6), and the conversion to the cc cross section from the sum of charm plus 
anticharm cross sections (0.5) [23]. We compare our total charm cross section 
to other experiments and to the NLO predictions as a function of pion-beam 
energy in Fig. 7. All experimental results are obtained by multiplying the 
D +D cross section by 1.7. The rise of the charm production cross section 
with energy is modeled reasonably well by the NLO theory, although the 
absolute value at any point depends greatly on the input parameters to the 
theory. 

In this paper we have presented the total forward cross section and differential 
cross sections versus Xp and pip for D mesons from Fermilab experiment E791 
data. This analysis represents the first measurement of the D cross section 
for a 500 GeV/c pion beam. The high statistics allows us to clearly observe 
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Fig. 6. The D +D p\ differential cross section compared to various theoretical 
predictions described in the text. The curves are normalized to obtain the best fit 
to the data in each case. The top plot shows the range 0<p^<4 (GeV/c) 2 while the 
bottom plot shows the full range, 0<p^<18 (GeV/c) 2 . Error bars do not include a 
1 15 % normalization uncertainty. 

a turnover point greater than zero (x c = 0.0131±0.0038) in the Feynman-rr 
distribution, providing evidence for a harder gluon distribution in the pion 
than in the nucleon. 

We have compared our differential cross section results to predictions from the 
next-to-leading order calculation by Mangano, Nason, and Ridolfi [12] and to 
the Monte Carlo event generator Pythia by T. Sjostrand et al. [9]. With suit- 
able choices for the intrinsic kt of the partons and the Peterson fragmentation 
function parameter, the NLO D meson calculation provides a good match to 
the py spectra and a fair match to the xp distribution. The string fragmenta- 
tion scheme in Pythia softens the original charm quark p^ distribution too 
much, and hardens the Xp spectra too much in both directions. However, the 
D result does predict the flattening of the xp cross section at high xp. The 
many adjustable parameters in the theoretical models allow one to obtain dis- 
tributions which are quite consistent with these data. Unfortunately, a given 



12 



Table 2 

Comparison of xf shape parameters to recent high-statistics pion-beam charm pro- 
duction experiments and to theory. E791 results are for D mesons; E769 (WA92) 
results are from a combined sample of D , D + , and D s (D and D + ) mesons. The 
theoretical results are obtained using the default parameters (described in the text). 
The fit range for n (Eq. 2) is in the table. The fit range for n' and x c (Eq. 3) is 
-0.125<x F <0.50. 



Experi- 


Energy 


xp Range 


n 


n! 


X c 


ment 


(GeV) 










E791 


500 


0.05-0.5 


4.61±0.19 


4.68±0.21 


.0131±.0038 


WA92[20] 


350 


0.0-0.8 


4.27±0.11 






E769[21] 


250 


0.0-0.8 


4.03±0.18 






MNR NLO c 


500 


0.05-0.5 


4.68 


5.06 


0.0231 


MNR NLO D 


500 


0.05-0.5 


5.53 


6.00 


0.0237 


Pythia c 


500 


0.05-0.5 


5.01 


5.12 


0.0115 


Pythia D° 


500 


0.05-0.5 


3.62 


3.66 


0.0041 



Table 3 

Comparison of p\ shape parameters to recent high-statistics pion-beam charm 
production experiments and to theory. Data samples are as described in Ta- 
ble 2. The fit range for b (Eq. 4) is 0<p^<4 (GeV/c) 2 except for WA92 which 
is 0<p?p<7 (GeV/c) 2 . The functions used to extract b' (Eq. 5) and a, (3 (Eq. 6) are 
fit in the range p\>\ (GeV/c) 2 and p^<18 (GeV/c) 2 , respectively. 



Experi- 


Energy 


b 


y 


a 


P 


ment 


(GeV) 


(GeV/c) 2 


(GeV/c)- 2 


(GeV/c)- 1 




E791 


500 


0.83±0.02 


2.41±0.03 


2.36±0.23 


5.94±0.39 


WA92[20] 


350 


0.89±0.02 








E769[21] 


250 


1.08±0.05 


2.74±0.09 


1.4±0.3 


5.0±0.6 


MNR NLO c 


500 


0.57 


1.88 


6.20 


8.68 


MNR NLO D 


500 


0.94 


2.32 


1.99 


5.30 


Pythia c 


500 


0.77 


2.09 


2.32 


5.14 


Pythia D° 


500 


1.06 


2.58 


1.55 


5.07 



set of parameters is neither unique, nor does it necessarily provide a good 
match to other data. In conjunction with other charm production results from 
this and other recent high-statistics experiments, however, it may be possi- 
ble to find a unique set of parameters. These results come from experiments 
with a variety of beam energies and types, and include measurements of dif- 
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Fig. 7. Theoretical and experimental results for the ir-N — ► cc cross section versus 
energy. The theoretical curves are obtained from MNR NLO predictions [23]. The 
three bands (solid, dashed, and dotted) correspond to three different charm quark 
masses (1.5, 1.2, and 1.8 GeV/c 2 ). The variation within the bands comes only from 
varying the renormalization scale. The factorization scale is kept fixed. All experi- 
mental data points are 1.7 times the single inclusive D +D xf>0 cross sections. 
The E653, NA27, WA92, NA32 (200 GeV), NA32 (230 GeV), and E769 data were 
obtained from references [24], [25], [20], [26], [27], and [28], respectively. 

ferential cross sections [20,28,31], production asymmetries [20,28-31,33], and 
correlations between two charm particles in the same event [3,35,36]. 

Unlike the uncertainties in the theoretical calculations of the differential cross 
sections, the uncertainties in the theoretical calculation of the total cross 
section come mostly from the perturbative calculation. The relatively large 
uncertainties are due to the low mass of the charm quark, which results in 
a large (unknown) contribution from higher-order terms. The total forward 
D°+D° cross section measured by E791 is a(D°+D° ; x F > 0) = 15.4 +£§ 
/ibarns/nucleon, assuming a linear atomic number dependence. The cross sec- 
tion is consistent with the MNR NLO prediction. 
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